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Abstract 

Background: Gaze processing deficits are a seminal, early, and enduring behavioral deficit in autism spectrum 
disorder (ASD); however, a comprehensive characterization of the neural processes mediating abnormal gaze 
processing in ASD has yet to be conducted. 

Methods: This study investigated whole-brain patterns of neural synchrony during passive viewing of direct and 
averted eye gaze in ASD adolescents and young adults {M Age = 16.6) compared to neurotypicals (NT) (M Age = 17.5) 
while undergoing magnetoencephalography. Coherence between each pair of 54 brain regions within each of three 
frequency bands (low frequency (0 to 15 Hz), beta (15 to 30 Hz), and low gamma (30 to 45 Hz)) was calculated. 

Results: Significantly higher coherence and synchronization in posterior brain regions (temporo-parietal-occipital) 
across all frequencies was evident in ASD, particularly within the low 0 to 15 Hz frequency range. Higher coherence in 
fronto-temporo-parietal regions was noted in NT. A significantly higher number of low frequency cross-hemispheric 
synchronous connections and a near absence of right intra-hemispheric coherence in the beta frequency band were 
noted in ASD. Significantly higher low frequency coherent activity in bilateral temporo-parieto-occipital cortical regions 
and higher gamma band coherence in right temporo-parieto-occipital brain regions during averted gaze was related 
to more severe symptomology as reported on the Autism Diagnostic Interview-Revised (ADI-R). 

Conclusions: The preliminary results suggest a pattern of aberrant connectivity that includes higher low frequency 
synchronization in posterior cortical regions, lack of long-range right hemispheric beta and gamma coherence, and 
decreased coherence in fronto-temporo-parietal regions necessary for orienting to shifts in eye gaze in ASD; a critical 
behavior essential for social communication. 
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Background 

Orienting to eye gaze is a vital skill present from birth and 
underlies effective non-verbal communication and social 
interaction [1] Failure to detect and/or respond in a 
typical manner to information conveyed by eye gaze is 
possibly the most important early hallmark of autism 
spectrum disorders (ASD) [2]. During social interaction, a 
persons eyes convey information about their direction of 
attention, emotion, and mental state [3]. Early responsive- 
ness to eye gaze cues develops quickly and evolves into 
joint attention, a complex form of non-verbal communica- 
tion that occurs when one individual follows another indi- 
viduals eye gaze or gesture to an object or third individual 
[1]. Joint attention is essential to the development of both 
social and language functioning [4-7], yet this pivotal skill 
fails to develop appropriately in infants with ASD [8,9]. 
Although research suggests that reflexive responses to eye 
gaze cues may be intact in individuals with ASD [10,11], 
studies indicate that orienting voluntarily to eye gaze is 
impaired and those with ASD may not demonstrate 
preferential sensitivity to eye gaze as a social cue [12-15]. 
Instead, it appears that for these individuals, others' eyes 
may merely serve as a spatial cue, much like a directional 
arrow [16]. 

Altered event-related potentials in response to shifts in 
eye gaze in 6- to 10-month-olds have been found to be 
associated with confirmed ASD diagnosis at 36 months 
of age, suggesting that abnormal response to eye gaze 
constitutes an important social cognitive endophenotype 
of ASD [17]. More specifically, in a recent investigation 
exploring neural responses to eye gaze cues in infants at 
risk for ASD, Elsabbagh and colleagues [17] reported 
that typically developing infants and infants at risk for 
but not diagnosed with ASD at 36 months of age dem- 
onstrated a higher P400 amplitude for gaze shifts away 
versus toward the infant, while those with a confirmed 
diagnosis at 36 months of age did not differentiate the 
conditions. The failure to appropriately orient to eye 
gaze in ASD likely has cascading effects on a wide range 
of developmental tasks including social and communica- 
tive functions [18]. However, while brain activation during 
tasks requiring eye gaze processing has been consistently 
found to be abnormal in ASD [19-23], the literature to 
date has not yet provided a coherent understanding of the 
neural processes underlying abnormal eye gaze processing 
in ASD. 

This study investigated whole-brain patterns of coher- 
ence during viewing of direct and averted eye gaze in ASD. 
Coherence, a measure of synchronization between active 
cortical network sources within a given frequency band 
[24], is used to measure short- and long-range connectivity. 
Theoretical models of ASD such as those proposed by 
Brock and colleagues [25,26] and Belmonte et al [27] 
propose that abnormally elevated levels of high-frequency 



neural activity and over-connectivity within localized brain 
regions causes impaired discrimination of signal from 
noise. This in turn causes impaired connectivity between 
distal cortical regions, leading to reduced activation of cor- 
tical regions involved in higher order processing compared 
to neurotypicals (NT). Recent research has supported these 
models, confirming that individuals with ASD consistently 
show abnormal patterns of connectivity between brain 
regions, with the most consistent finding being a lower 
degree of connectivity between frontal and posterior 
regions compared to NT, both at rest and while performing 
cognitive tasks [28]. 

Recent reports suggest that disruptions in synchronous 
neural oscillatory activity are a primary cellular mechanism 
of impairment in ASD [29], particularly in the gamma 
frequency band, as identified by resting-state EEG and 
magnetoencephalography (MEG) [15-17]. Gamma band 
synchrony has been found to be associated with percep- 
tual binding at early levels of sensory processing, attention 
[30], and working and long-term memory [31], and is pur- 
ported to be involved in top-down modulation of sensory 
signals and large-scale integration of distributed neural 
networks [32], 

Despite recent interest in neural connectivity in ASD, 
very few MEG studies have examined functional connect- 
ivity during face processing, and no studies have examined 
connectivity during eye gaze processing. In a recent MEG 
investigation in children and adults with ASD, the gamma 
frequency band response in right lateral occipital areas 
was largely absent compared to typically developing 
participants while viewing emotions on faces [33]. Khan 
and colleagues also described abnormalities in gamma os- 
cillatory activity during face processing in ASD compared 
to typically developing adolescents and young adults [33] . 
The authors focused on nesting oscillations in which the 
amplitude or phase of a lower frequency band (for 
example, alpha) modulates the phase or amplitude of a 
higher frequency band, referred to as phase amplitude 
coupling (PAC). Group differences in alpha-gamma PAC 
were noted but in the absence of reductions in alpha or 
gamma power, suggesting these differences were driven by 
variation in the timing of gamma generation. In addition, 
statistically significant differences in coherence in the 
alpha band in the inferior frontal gyrus, anterior cingulate, 
and left precuneus were noted, consistent with reductions 
in long-range connectivity. Reductions in local functional 
connectivity within the fusiform face area (FFA) were also 
reported. 

The central role of eye gaze processing impairment in 
ASD makes eye gaze a good candidate for elucidating 
aberrant patterns of connectivity and as a potential 
biomarker. We utilized MEG to record patterns of neural 
activation during direct and averted eye gaze processing in 
ASD and NT, enabling calculation of coherence between 
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circumscribed cortical regions within discrete frequency 
bands. Specifically, we calculated coherence within the beta 
and gamma frequency bands, which are currently impli- 
cated in long- and short-range transmission of information 
involved in high-level cognitive processing, respectively 
[34]. As an exploratory aim, we also examined a lower' 
frequency band, collapsing delta, theta, and alpha bands 
for comparison with the alert working brain, beta, and 
gamma frequencies. 

Methods 

Participants 

Eighteen participants completed the study; 10 participants 
with ASD (M Age = 16.6); M IQ = 112) and eight neurotypical 
controls (M Age = 17; M IQ =116). The groups did not differ 
significantly in age {11(16) = 1.79, P = 0.76), gender (/ 2 = 
0.11), or Full Scale IQ, with both groups generally 



performing in the Average to Above Average range on 
the Wechsler Abbreviated Scale of Intelligence (WASI; 
£/(16) = 40, P = 0.74) [35]. See Table 1 for a review of 
the demographic variables. There was equal age distri- 
bution between the genders of the two groups. To en- 
sure that potential developmental differences between 
the groups did not account for our effect, significant re- 
lationships between our brain region coherence values 
with age were inspected. Age was related to coherence 
values in the beta frequency band only in the right 
inferior frontal to right superior temporal pathway in 
the NT controls (r = 0.52, P = 0.02). Five of the 10 ASD 
participants were on psychotropic medications. Two 
participants were prescribed a single medication, an 
antidepressant (Prozac) or antipsychotic (Risperdal). 
Two participants were prescribed an antidepressant 
and psychostimulant (Prozac or Zoloft and Concerta), 



Table 1 Demographic characteristics of the participants 



ASD subjects 


Gender 


Age 


Dominance 


FSIQ standard 
score 


Vocabulary 
T-score 


Matrix reasoning 
T-score 


ADI-R 
social 


ADI-R 

communication 


ADIR 

repetitive 

behaviors 


1 


M 


19.00 


R 


131.0 


77.00 


57.00 


14.00 


16.00 


9.00 


2 


F 


13.00 


L 


96.0 


35.00 


61.00 


23.00 


15.00 


4.00 


3 


M 


13.00 


R 


93.0 


39.00 


53.00 


27.00 


24.00 


6.00 


4 


M 


15.00 


R 


137.0 


67.00 


73.00 


20.00 


14.00 


8.00 


5 


M 


19.00 


R 


144.0 


77.00 


67.00 


17.00 


9.00 


3.00 


6 


F 


16.00 


R 


119.0 


62.00 


60.00 


20.00 


11.00 


5.00 


7 


M 


24.00 


R 




Discontinued 


66.00 


24.00 


18.00 


3.00 


8 


F 


16.00 


R 


79.0 


44.00 


28.00 


20.00 


18.00 


6.00 


9 


M 


16.00 


R 


118.0 


61.00 


60.00 


13.00 


10.00 


7.00 


10 


M 


15.00 


R 


92.0 


64.00 


28.00 


17.00 


16.00 


8.00 


Mean 




16.60 




112.1 


58.44 


55.30 


19.50 


15.10 


5.90 


SD 




3.31 




22.9 


15.61 


15.42 


4.40 


4.46 


2.13 


Range 




11-24 




79-144 


35-77 


28-73 








Control subjects 


1 


M 


21.00 


R 


122.0 


64.00 


61.00 








2 


F 


13.00 


R 


120.0 


66.00 


57.00 








3 


M 


13.00 


R 


122.0 


59.00 


66.00 








4 


M 


15.00 


R 


118.0 


60.00 


60.00 








5 


M 


19.00 


R 


117.0 


58.00 


60.00 








6 


F 


19.00 


L 


119.0 


59.00 


63.00 








7 


M 


18.00 


R 


100.0 


44.00 


57.00 








8 


F 


18.00 


R 


117.0 


60.00 


57.00 








Mean 




17.00 




116.9 


58.75 


60.13 








SD 




2.98 




7.1 


6.56 


3.23 








Range 




13-21 




100-122 


44-66 


57-66 








Mann-Whitney U P value 


0.11 


0.76 




0.70 


0.54 


0.90 









ADI-R = Autism Diagnostic Interview-Revised; ASD = autism spectrum disorder; FSIQ = Full Scale Intelligence Quotient. 
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and one participant was prescribed three medications, an 
antidepressant, anxioltyic, and mood stabilizer (Cymbalta, 
Buspar, and Depakote). 

Individuals were recruited from and underwent MEG 
procedures at Henry Ford Hospital (HFH). Participants 
were diagnosed with Pervasive Developmental Disorder 
(PDD) (recently revised to ASD) based on the Diagnostic 
and Statistical Manual of Mental Disorders-Fourth 
Edition-Text Revision (DSM-IV-TR) [2] diagnostic 
criteria. Diagnoses were confirmed with the Autism 
Diagnostic Interview-Revised (ADI-R) [36]. Means and 
standard deviations for the ADI-R domain scores include: 

^social (SD) = 19.25 (4.59); ^Communication (SD) = 15.50 

(4.34); M Rep etitive (SD) = 6.63 (1.69). Inclusion criteria in- 
cluded at least low average intelligence (>80 Full Scale IQ 
scores on the WASI). Exclusionary criteria included any 
known history of head injury with unconsciousness, epi- 
lepsy, affective or anxiety disorders. No ASD participant 
had a history of a genetic disorder. NT participants had no 
history of developmental delay, learning disorder, or ASD 
in a first-degree relative. All APA Ethical Guidelines were 
followed and Institutional Review Board approval was 
obtained from all institutions participating in this study. 

MEG data acquisition and preprocessing 

Cortical activity was recorded using a 148 channel whole 
head MEG system (4D Neuroimaging, Magnes WH2500) 
with magnetometer type sensors. During acquisition, the 
data were band-pass filtered 0.1 to 100 Hz and digitally 
sampled at 508.63 Hz. The timing of stimuli was recorded 
as pulse codes (representing the type of stimulus) on a 
trigger channel simultaneously collected with the MEG 
data. In postprocessing, noise artifacts due to heart and 
body movement were eliminated using an independent 
component analysis (ICA). Singular valued decomposition 
was used to remove any other artifacts in the data such as 
mouth movements if needed. Regarding movement 
artifact, runs are repeated if the coil on head positions ex- 
ceeds 0.5 cm, although this did not occur during acquisi- 
tion. As such, we combined the runs before source 
reconstruction. Regarding noise reduction, 4D Neuroim- 
aging incorporates a set of reference sensors that are used 
to sample the environmental magnetic fields and create a 
file that has a set of weights. These weights are subtracted 
from the data during data collection. Data were band-pass 
filtered from 1 to 50 Hz. The locations of events on the 
trigger and response channels were used to select 2-s 
epochs of MEG data to examine average evoked responses 
for the stimuli requiring a conditional button press (that 
is, asterisk, face, or words). All trials within each condition 
were averaged to determine the evoked response. All 
epochs had a baseline of 500 ms before stimuli onset and 
1,500 ms of data after stimuli onset. 



Gaze cueing paradigm 

The gaze cueing paradigm used was an adaptation of a 
paradigm used by Pelphrey and colleagues [37]. MEG field 
responses to gaze cues were collected for two 14-min tri- 
als in which five task conditions were administered: direct 
gaze, averted gaze, and gaze cueing to peripheral stimuli 
(asterisk, word, or face). In the direct and averted gaze 
conditions, participants passively viewed a central charac- 
ter. There were 30 trials in the direct gaze condition and 
30 each (left and right) in the averted gaze condition. In 
each of the three additional gaze cueing conditions, the 
central character engaged in a random gaze shift toward 
the right or left for 1 s. A target (asterisk, word, or face) 
then appeared at either the right or the left of the subject 
for 3 s. See Figure 1 for an example of one of the target 
stimuli. The next trial began with the character returning 
to a direct gaze for 2 s with no stimuli in the periphery. 
The location of the target stimulus was either congruent 
or incongruent with the direction of the character s gaze. 
Sixty targets were presented in each gaze-cueing condi- 
tion, including 30 congruent and 30 incongruent trials. A 
conditional button press during gaze cues to the periph- 
eral stimuli conditions ensured engagement during the 
passive conditions. For the purposes of this analysis, only 
passive viewing of direct and averted gaze was considered. 
Results from the evoked response data requiring active re- 
sponses to the conditional stimuli (asterisk, faces, words) 
have been previously published [38]. 

MEG coherence analysis 

Synchronization of neuronal activity was quantified by 
calculating coherence between cortical sites from MEG 
imaged brain activations [39,40]. A model of the cortical 
brain surface was created from an age- and gender- 
appropriate standard MRI. The MRI was segmented and 
the brain surface was represented by a cortical model of 
approximately 4,000 dipoles each having an x, y, and z 
orientation at each site. Sites were distributed to represent 
the same volume of cortical gray matter. This model was 
then morphed to fit the digitized head shape collected 
during the MEG acquisition. To calculate coherence [40], 
the MEG data were first divided into 80 segments each 
containing 7.5-s segments of data and cortical activity in 
each segment was imaged on to the MRI using the MR- 
FOCUSS imaging technique [41]. Using the time sequence 
of imaged activity, coherence between active cortical 
model sites was calculated for each data segment and then 
averaged for the completed study. In addition, for each 
cortical model site, connectivity was quantified by a histo- 
gram of the number of sites to which the site had the 
same level of coherence. Statistical analysis of cortical co- 
herence levels (0 to 1) were used to quantify differences in 
network connectivity between groups. Changes in coher- 
ence and connectivity between brain regions implicated as 
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Figure 1 Gaze paradigm. Analyses included passive conditions corresponding to slides with direct and averted gaze. In this example, the target 
slide would require a button press as the character's gaze is looking at the asterisk. This was how attention to eye gaze was ensured. 



having deviant electrophysiological activity in the ASD 
brain were quantified and subjected to further statistical 
analysis. 

Power spectra for activity at all active sites were also cal- 
culated and used to quantify differences in low frequency, 
beta, and gamma power. For the 'lower' frequency band, 
delta, theta, and alpha bands were collapsed for compari- 
son with the alert working brain, specifically beta and 



gamma frequencies. For additional details of MEG coher- 
ence imaging, see our publication [40]. A region-of- 
interest (ROI) tool implemented in MEG Tools was used 
to identify 54 regions in the brain (27 in each hemisphere). 
See Figure 2 for a list of the brain regions. MEG Tools 
uses a non-linear volumetric transformation of the pa- 
tients brain to transform MEG coordinates to standard 
Talairach [42] or MM [43] coordinates [44]. This enables 




Figure 2 Brain regions used to develop 1,431 brain region pathways. Total brain regions equal 54: 27 right and 27 left hemisphere. 

(A) Midsagittal, (B) sagittal, (C) lateral (postmortem), and (D) coronal MR-T1 weighted images revealing 25 of the 27 brain regions. Gyrus rectus 

and middle orbitofrontal were unable to be visualized given that all x, y, and z locations cannot be viewed concurrently. 
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the ROI tool to access an atlas of Brodmann s area identi- 
fiers and an atlas of cortical structures [43] . 

Group difference testing 

For each frequency band (low, beta, and gamma) within 
each condition (direct and averted gaze), a £-test was con- 
ducted to assess group difference in average coherence 
values for each pair of brain regions (N = 1,431). A P value 
was produced for each region pair. The false discovery 
rate (FDR) was used to adjust for multiple testing. Because 
of the large number of tests being performed simultan- 
eously, using a significance level of alpha = 0.05 without 
adjusting for multiple testing would lead to a large num- 
ber of false positive results. Bonferroni adjustments for 
multiple comparisons aim to control the Family Wise 
Error Rate. If a Bonferroni correction were applied to 
every test, there would be only a 5% chance of at least one 
false positive in all of the N = 1,431 tests. Bonferroni cor- 
rections require the P value to be less than 0.05/N, where 
N is the number of tests. With N = 1,431, this criterion 
becomes especially stringent, and many true differences 
may be missed (false negatives). In the era of large scale 
testing, a less conservative approach to adjusting for mul- 
tiple testing has been developed. The FDR is the propor- 
tion of tests declared significant that are actually different 
only due to chance (or the proportion of significant tests 
that are false positives). The FDR is a widely accepted, less 
conservative approach to adjusting for multiple testing in 
large scale problems. The Benjamini-Hochberg algorithm 
[45] was used to control the FDR at 0.10. From each £-test, 
a z-score was computed according to the method of Efron 
[46] to summarize the difference in coherence values 
between ASD and NT. Positive ^-scores indicate higher 
coherence in the ASD group. A series of chi-squares were 
computed to determine if the number of intra-hemispheric 
and inter-hemispheric cortical differences within the low (0 
to 15 Hz), beta (15 to 30 Hz), and low gamma (30 to 
45 Hz) frequency bands were statistically different between 
the groups. 

To examine relationships between ASD clinical symp- 
toms as reported on the ADI-R and neural oscillatory 
activity (coherence), a series of Kendall Tau correlation 
coefficients were computed. 

Results 

Behavioral data 

There were no significant between group differences noted 
in error rates in responding to the conditional button press 
during the intervening task condition with respect to ac- 
curacy (£(16) = 0.70 P = 0.51), suggesting that both groups 
were equally engaged in the task. Reaction times were not 
statistically different (£(16) = -0.11 P = 0.92). See Table 2 for 
the accuracy (Total Correct Responses) and reaction time 
values for correct responses for the NT and ASD groups. 



MEG analysis: coherence imaging of connectivity 
Direct gaze condition 

During the direct gaze condition, 91 of the 1,431 pathways 
were found to be significantly different between the 
groups. In NT, higher coherence was observed between 
frontal, temporal, and parietal regions. Higher coherence 
was particularly evident between bilateral frontal (middle, 
inferior, and orbitofrontal) gyri and right superior tem- 
poral, pre- and postcentral gyri. In ASD participants, 
higher coherence was noted between left occipito-parietal 
(angular, middle, and superior occipital gyri) and bilateral 
occipito-parietal regions (inferior, middle, superior occipi- 
tal gyri, and supramarginal regions). See Table 3 for brain 
region pairs with statistically significant between group 
differences in coherence for direct gaze collapsed across 
frequency bands. Only the top 10 regions for each group 
are presented to simplify the results. No significant group 
differences were found within the separate frequency 
bands during the direct gaze condition. 

Averted gaze condition 

During the averted gaze condition, 390 of the 1,431 path- 
ways were found to be significantly different between the 
groups. Consistent with the direct gaze findings, NT dem- 
onstrated significantly higher coherent activity across all 
frequencies in fronto-temporo-parietal regions, consistent 
with known neuroanatomical substrates critical for 
responding to shifts in eye gaze (see Figure 3). That is, sig- 
nificantly higher coherence was noted between bilateral 
frontal (inferior, middle, superior, orbitofrontal gyri) and 
right frontal (inferior, middle, superior, and precentral 
gyri), superior temporal, and parietal (postcentral gyrus) 
regions. ASD participants displayed higher coherence be- 
tween left parieto-occipital (angular, inferior, and middle 
occipital) and bilateral temporo-parieto-occipital regions 
(inferior, middle, superior temporal, occipital, angular gyri) 
(see Figure 3). 

In contrast to direct gaze, statistically significant be- 
tween group differences were noted within each frequency 
during the averted gaze condition. Of the 390 pathways, 
significant differences in coherence were found within 
specific frequencies in 233 of the pathways; 127 in the low 
frequency band, 37 in the beta frequency band, and 69 in 
the low gamma frequency band. Table 4 presents the brain 
region pairs with statistically significant between group 
differences in coherence for averted gaze in the low (0 to 
15 Hz), beta (15 to 30 Hz), and gamma (30 to 45 Hz) 
frequency bands. 

In the low frequency band, ASD participants displayed 
higher coherent activity between left parieto-occipital 
regions and right temporo-parieto-occipital regions and 
significantly lower coherence between bilateral frontal and 
right fronto-temporo-parietal regions. In the beta band, 
ASD participants demonstrated higher coherence between 
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Table 2 Average accuracy (total correct responses) and reaction time values for correct responses for each ASD and 
neurotypical participant 



Total correct responses 
ASD 



Neurotypical 



Reaction time for correct responses 
ASD 



Neurotypical 



47 
43 
46 
37 
45 
38 
53 
45 
45 
41 

M (SD) 44.00 (4.62) 



45 
46 
47 
44 
45 
45 
45 
38 



M (SD) 44.38 (2.72) 



673 

1,044 

818 

1,096 

683 

936 

935 

761 

557 

1,136 

M (SD) 863.90 (196.26) 



915 
916 
887 
717 
718 
855 
862 
929 



M (SD) 849.88 (88.15) 



M = Mean; SD = Standard Deviation. 
Values reflect averages over the two runs. 



Table 3 Brain region pairs with statistically significant between group differences in coherence for direct gaze 
collapsed across frequency bands 

Brain region pairs (pathway) for direct gaze z-score P value 

Regions with higher coherence in NT 

R. middle frontal gyrus and R. superior temporal gyrus -4.69 <0.0001 

L. lateral orbitofrontal gyrus and R. precentral gyrus -4.44 <0.0001 

L. lateral orbitofrontal gyrus and R. postcentral gyrus -4.41 <0.0001 

R. inferior frontal gyrus and R. middle frontal gyrus -4.33 <0.0001 

L. superior frontal gyrus and R. superior temporal gyrus -4.32 <0.0001 

L. inferior frontal gyrus and R. postcentral gyrus -4.26 <0.0001 

L. inferior frontal gyrus and R. precentral gyrus -4.14 <0.0001 

R. middle frontal gyrus and R. postcentral gyrus -4.19 <0.0001 

L. middle frontal gyrus and R. precentral gyrus -4.1 1 0.0001 

L. middle frontal gyrus and R. postcentral gyrus -4.12 0.0001 
Regions with higher coherence in ASD 

L. angular gyrus and L. middle occipital gyrus 5.10 <0.0001 

L. angular gyrus and R. middle occipital gyrus 4.77 <0.0001 

L. angular gyrus and L. superior occipital gyrus 4.09 <0.0001 

L. angular gyrus and R. cuneus 4.02 0.0001 

L. angular gyrus and L. inferior occipital gyrus 4.02 0.0001 

L. postcentral gyrus and L. superior occipital gyrus 3.86 0.0001 

L. angular gyrus and R. inferior occipital gyrus 3.85 0.0001 

L. superior occipital gyrus and L. supramarginal gyrus 3.85 0.0001 

L. middle occipital gyrus and R. middle occipital gyrus 3.84 0.0001 

L. angular gyrus and L. supramarginal gyrus 3.71 0.0002 



Only top 10 for each group of the 91 total are reported to simplify results. 
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Occi p i tal- Fron la I 



Figure 3 Subset of regions of network of greatest average coherence in a representative ASD and NT participant. Red reflects sending of 
information while green reflects receiving of information flow in areas of high average coherence. In a male adolescent with ASD (A), information in 
posterior cortical areas (parietal cortex) is highly active in sending information when processing eye gaze cues while minimal frontal activity is noted. 
In a NT adolescent male (B), the left inferior and middle frontal gyri are highly active in sending information when processing gaze cues. 



left parietooccipital regions and bilateral temporo-occipital 
and left parietal regions. In the gamma band, ASD partici- 
pants showed higher coherence between bilateral temporo- 
parietooccipital regions as well as bilateral parietal and 
orbitofrontal regions. In both the beta and gamma fre- 
quency bands, ASD participants showed lower coherence 
between bilateral frontal, fronto-temporal, and temporo- 
parietal regions compared to NT (See Figure 4). In contrast, 
NT subjects displayed significantly higher coherence be- 
tween bilateral frontal, fronto-temporal, and fronto-parietal 
regions across all frequency bands. Due to the essentially 
non-overlapping distributions, formal chi-square or log 
linear analyses could not be computed. The number of 
region-to-region cortical connection counts within the low 
(0 to 15 Hz), beta (15 to 30 Hz), and low gamma (30 to 
45 Hz) frequency bands can be visually inspected in 
Figure 4. 

Regarding intra- and inter-hemispheric differences in 
coherence, there was a significant association between 
group membership and the number of coherent right 
intra-hemispheric connections X 2 (2) = 8.34, P <0.01. Es- 
sentially no right intra-hemispheric coherent connections 
were noted in the 15 to 30 Hz range during passive view- 
ing of averted gaze in ASD. No association between group 
membership and the number of coherent left intra- 
hemispheric connections was noted X 2 (2) = 2.57, P = 027 
There was also a significant association between group 
membership and the number of inter-hemispheric coher- 
ent connections X 2 (2) = 10.01, P <0.007. More specifically, 
the ASD participants had a 6.5 times greater number of 
low frequency inter-hemispheric coherent connections 
relative to coherent beta and gamma frequency connec- 
tions than NT while viewing averted gaze. See Figure 5 to 
examine differences in intra- and inter-hemispheric cor- 
tical connections within the low, beta, and low gamma 



frequency bands in regions of either higher coherence in 
NT or ASD participants during averted gaze. 

Relationship between oscillatory activity and clinical 
symptomatology 

ADI-R scores were significantly correlated with coher- 
ence values in a predictable pattern among the ASD 
group under both the direct and averted gaze conditions 
(see Table 5). During direct gaze and collapsed over all 
frequencies, significantly higher coherence between left 
parietal regions (that is, angular and supramarginal gyri) 
was related to a greater number of ASD symptoms as re- 
ported on the ADI-R (<0.005) while higher coherence 
between bilateral frontal cortical regions was related to 
fewer ASD symptoms. That is, between left inferior 
frontal/lateral orbitofrontal and right precentral gyri 
(<0.01 and <0.005, respectively). 

A similar pattern of relationships was noted during 
the averted gaze condition in the low frequency band. 
Significantly higher coherence between unilateral and bi- 
lateral parietal regions (P <0.005), between left parietal 
and left occipital cortices (P <0.005) as well as between 
left parietal and bilateral temporal regions (P <0.005), 
and between left temporal and bilateral occipital regions 
(R P <0.005; L P <0.01) was related to a higher scores on 
the ADI-R and more severe ASD symptoms. In con- 
trast, higher coherence between unilateral right frontal 
{P <0.001) and bilateral frontal regions (P <0.005) as well 
as between right frontal and right parietal (P <0.005) and 
right frontal and right occipital (P <0.005) was related to 
fewer ASD symptoms. 

Primarily negative correlations were obtained during the 
averted gaze condition in the beta frequency band. Higher 
beta frequency band coherence between bilateral frontal 
regions (P <0.005), between unilateral right frontal regions 
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Table 4 Brain region pairs with statistically significant between group differences in coherence for averted gaze in the 
(A) low (0 to 15 Hz), (B) beta (15 to 30 Hz), and (C) gamma (30 to 45 Hz) frequency bands 

(A) Brain Region Pairs (Pathway) for 0 to 15 Hz z-score P value 

Regions with higher coherence in NT 



R. inferior frontal gyrus and R. middle frontal gyrus 


-5.23 


<0.0001 


R. inferior frontal gyrus and R. superior frontal gyrus 


-4.96 


<0.0001 


L. lateral orbitofrontal gyrus and R. postcentral gyrus 


-4.2 


<0.0001 


R. superior frontal gyrus and R. superior temporal gyrus 


-4.2 


<0.0001 


R. postcentral gyrus and R. superior frontal gyrus 


-4.07 


0.00010 


L. precentral gyrus and R. superior frontal gyrus 


-4.02 


0.00020 


L. superior frontal gyrus and R. postcentral gyrus 


-3.65 


0.00030 


L. inferior frontal gyrus and R. postcentral gyrus 


-3.65 


0.00040 


L. lateral orbitofrontal gyrus and R. middle frontal gyrus 


-3.61 


0.00030 


L. inferior frontal gyrus and R. inferior frontal gyrus 


-3.61 


0.00030 


L. superior temporal gyrus and R. middle frontal gyrus 


-3.55 


0.00040 


L. precentral gyrus and R. middle frontal gyrus 


-3.52 


0.00080 


R. middle frontal gyrus and R. superior temporal gyrus 


-3.51 


0.00060 


L. lateral orbitofrontal gyrus and R. supramarginal gyrus 


-3.49 


0.00060 


R. middle frontal gyrus and R. postcentral gyrus 


-3.49 


0.00080 


L. middle frontal gyrus and R. postcentral gyrus 


-3.42 


0.00100 


L. lateral orbitofrontal gyrus and R. inferior frontal gyrus 


-3.3 


0.00130 


L. middle frontal gyrus and R. inferior frontal gyrus 


-3.27 


0.00110 


L. middle frontal gyrus and R. precentral gyrus 


-3.27 


0.00150 


R. superior frontal gyrus and R. supramarginal gyrus 


-3.26 


0.00170 


L. lateral orbitofrontal gyrus and R. precentral gyrus 


-3.25 


0.00130 


L. superior temporal gyrus and R. superior frontal gyrus 


-3.24 


0.00120 


L. precentral gyrus and L. superior frontal gyrus 


-3.22 


0.00190 


R. precentral gyrus and R. superior frontal gyrus 


-3.2 


0.00150 


R. middle frontal gyrus and R. superior frontal gyrus 


-3.14 


0.00170 


L. lateral orbitofrontal gyrus and L. precentral gyrus 


-3.14 


0.00220 


L. superior frontal gyrus and R. inferior frontal gyrus 


-3.1 


0.00200 


L. inferior frontal gyrus and R. precentral gyrus 


-3.08 


0.00250 


R. angular gyrus and R. superior frontal gyrus 


-3.07 


0.00320 


L. inferior frontal gyrus and R. middle frontal gyrus 


-3.05 


0.00230 


L. postcentral gyrus and R. superior frontal gyrus 


-3.02 


0.00320 


L. inferior frontal gyrus and R. supramarginal gyrus 


-3.02 


0.00320 


L. superior frontal gyrus and R. precentral gyrus 


-3 


0.00270 


R. middle frontal gyrus and R. supramarginal gyrus 


-2.92 


0.00470 


L. superior parietal gyrus and R. superior frontal gyrus 


-2.88 


0.00490 


L. precuneus and R. superior frontal gyrus 


-2.83 


0.00730 


L. inferior frontal gyrus and L. precentral gyrus 


-2.82 


0.00560 


L. superior frontal gyrus and R. supramarginal gyrus 


-2.81 


0.00610 


L. lateral orbitofrontal gyrus and R. angular gyrus 


-2.8 


0.00570 


L. lateral orbitofrontal gyrus and L. superior parietal gyrus 


-2.78 


0.00650 


L. middle orbitofrontal gyrus and R. postcentral gyrus 


-2.78 


0.00700 


L. middle frontal gyrus and R. supramarginal gyrus 


-2.78 


0.00720 



Lajiness-O'Neill et al. Journal of Neurodevelopmental Disorders 2014, 6:15 
http://www.jneurodevdisorders.eom/content/6/1 /1 5 



Page 10 of 22 



Table 4 Brain region pairs with statistically significant between group differences in coherence for averted gaze in the 
(A) low (0 to 15 Hz), (B) beta (15 to 30 Hz), and (C) gamma (30 to 45 Hz) frequency bands (Continued) 


L. lateral orbitofrontal gyrus and L. precuneus 


-2.76 


0.00760 


L. superior temporal gyrus and R. inferior frontal gyrus 


-2.72 


0.00660 


L. middle frontal gyrus and L. precentral gyrus 


-2.71 


0.00880 


R. middle frontal gyrus and R. precentral gyrus 


-2.68 


0.00850 


L. superior frontal gyrus and R. superior temporal gyrus 


-2.66 


0.00780 


Regions with higher coherence in ASD 






L. angular gyrus and L. inferior occipital gyrus 


4.86 


<0.0001 


L. angular gyrus and R. lingual gyrus 


4.79 


<0.0001 


L. middle occipital gyrus and R. cuneus 


4.7 


<0.0001 


L. middle occipital gyrus and R. lingual gyrus 


4.51 


<0.0001 


L. angular gyrus and L. middle occipital gyrus 


4.47 


<0.0001 


L. angular gyrus and R. middle occipital gyrus 


4.46 


<0.0001 


L. angular gyrus and R. cuneus 


4.38 


<0.0001 


L. inferior occipital gyrus and R. lingual gyrus 


4.28 


<0.0001 


L. angular gyrus and R. superior occipital gyrus 


4.17 


<0.0001 


L. middle occipital gyrus and R. superior occipital gyrus 


4.16 


<0.0001 


L. middle occipital gyrus and R. middle occipital gyrus 


4.11 


<0.0001 


L. inferior occipital gyrus and R. cuneus 


4.01 


0.00010 


L. cuneus and R. lingual gyrus 


3.99 


0.00010 


L. middle occipital gyrus and R. superior parietal gyrus 


3.9 


0.00010 


R. lingual gyrus and R. middle occipital gyrus 


3.83 


0.00010 


L. inferior occipital gyrus and R. superior parietal gyrus 


3.74 


0.00020 


L. angular gyrus and R. inferior occipital gyrus 


3.72 


0.00020 


L. angular gyrus and L. cuneus 


3.68 


0.00020 


R. cuneus and R. middle occipital gyrus 


3.6 


0.00040 


L. middle occipital gyrus and R. inferior occipital gyrus 


3.59 


0.00040 


R. gyrus rectus and R. lingual gyrus 


3.57 


0.00060 


L. angular gyrus and L. lingual gyrus 


3.55 


0.00040 


L. middle temporal gyrus and R. lingual gyrus 


3.54 


0.00040 


L. angular gyrus and L. fusiform gyrus 


3.54 


0.00060 


R. cuneus and R. lingual gyrus 


3.51 


0.00060 


R. middle occipital gyrus and R. superior occipital gyrus 


3.47 


0.00050 


L. middle occipital gyrus and R. angular gyrus 


3.47 


0.00050 


L. middle occipital gyrus and R. middle temporal gyrus 


3.47 


0.00050 


L. superior temporal gyrus and R. cuneus 


3.45 


0.00070 


L. middle temporal gyrus and R. cuneus 


3.44 


0.00070 


L. lingual gyrus and R. lingual gyrus 


3.41 


0.00070 


L. inferior temporal gyrus and R. cuneus 


3.39 


0.00070 


L. inferior occipital gyrus and R. middle occipital gyrus 


3.34 


0.00090 


L. angular gyrus and L. inferior temporal gyrus 


3.34 


0.00110 


R. lingual gyrus and R. superior occipital gyrus 


3.32 


0.00090 


L. inferior occipital gyrus and R. superior occipital gyrus 


3.31 


0.00090 


L. postcentral gyrus and R. cuneus 


3.3 


0.00100 


L. fusiform gyrus and R. cuneus 


3.3 


0.00110 
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Table 4 Brain region pairs with statistically significant between group differences in coherence for averted gaze in the 
(A) low (0 to 15 Hz), (B) beta (15 to 30 Hz), and (C) gamma (30 to 45 Hz) frequency bands (Continued) 



L. angular gyrus and R. superior parietal gyrus 


3.26 


0.00120 


L. inferior temporal gyrus and R. lingual gyrus 


3.25 


0.00120 


L. inferior occipital gyrus and R. middle temporal gyrus 


3.23 


0.00140 


R. middle occipital gyrus and R. superior parietal gyrus 


3.21 


0.00130 


R. cuneus and R. superior occipital gyrus 


3.19 


0.00150 


L. cuneus and R. middle occipital gyrus 


3.19 


0.00160 


L. angular gyrus and L. middle temporal gyrus 


3.17 


0.00160 


L. superior temporal gyrus and R. lingual gyrus 


3.16 


0.00170 


L. angular gyrus and L. superior occipital gyrus 


3.12 


0.00190 


L. cuneus and R. cuneus 


3.1 


0.00230 


R. cuneus and R. gyrus rectus 


3.09 


0.00240 


L. cuneus and R. inferior occipital gyrus 


3.08 


0.00210 


L. angular gyrus and R. middle temporal gyrus 


3.04 


0.00250 


L. lingual gyrus and R. cuneus 


3.02 


0.00280 


L. cuneus and L. fusiform gyrus 


3.01 


0.00290 


L. inferior occipital gyrus and R. angular gyrus 


2.98 


0.00300 


L. cuneus and L. middle occipital gyrus 


2.98 


0.00380 


L. middle occipital gyrus and R. precuneus 


2.95 


0.00320 


L. postcentral gyrus and R. lingual gyrus 


2.94 


0.00340 


L. inferior occipital gyrus and L. middle occipital gyrus 


2.94 


0.00480 


L. postcentral gyrus and R. superior occipital gyrus 


2.88 


0.00460 


R. superior occipital gyrus and R. superior parietal gyrus 


2.87 


0.00410 


L. inferior occipital gyrus and R. inferior occipital gyrus 


2.86 


0.00450 


L. angular gyrus and R. angular gyrus 


2.85 


0.00440 


L. lingual gyrus and R. middle occipital gyrus 


2.83 


0.00480 


L. fusiform gyrus and R. lingual gyrus 


2.83 


0.00490 


L. middle occipital gyrus and R. lateral orbitofrontal gyrus 


2.81 


0.00500 


L. postcentral gyrus and L. superior occipital gyrus 


2.81 


0.00530 


R. lingual gyrus and R. superior parietal gyrus 


2.78 


0.00540 


L. cuneus and R. superior parietal gyrus 


2.78 


0.00560 


L. middle occipital gyrus and L. middle temporal gyrus 


2.76 


0.00750 


L. middle occipital gyrus and R. supramarginal gyrus 


2.75 


0.00630 


L. precentral gyrus and R. cuneus 


2.74 


0.00630 


R. lingual gyrus and R. middle temporal gyrus 


2.74 


0.00680 


R. inferior occipital gyrus and R. lingual gyrus 


2.73 


0.00650 


L. cuneus and R. superior occipital gyrus 


2.7 


0.00740 


L. angular gyrus and R. gyrus rectus 


2.69 


0.00740 


R. middle occipital gyrus and R. middle temporal gyrus 


2.69 


0.00740 


L. angular gyrus and L. superior temporal gyrus 


2.68 


0.00740 


L. cuneus and L. postcentral gyrus 


2.67 


0.00770 


L. angular gyrus and R. precuneus 


2.67 


0.00790 


L. middle occipital gyrus and R. superior temporal gyrus 


2.64 


0.00850 


(B) Brain region pairs (pathway) for 15 to 30 Hz 


z-score 


P value 



Regions with higher coherence in NT 
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Table 4 Brain region pairs with statistically significant between group differences in coherence for averted gaze in the 
(A) low (0 to 15 Hz), (B) beta (15 to 30 Hz), and (C) gamma (30 to 45 Hz) frequency bands (Continued) 



R. middle frontal gyrus and R. superior temporal gyrus 

R. inferior frontal gyrus and R. superior frontal gyrus 

L. lateral orbitofrontal gyrus and R. precentral gyrus 

R. inferior frontal gyrus and R. superior temporal gyrus 

L. superior frontal gyrus and R. precentral gyrus 

R. middle frontal gyrus and R. precentral gyrus 

R. middle frontal gyrus and R. postcentral gyrus 

L. superior frontal gyrus and R. superior temporal gyrus 

L. superior frontal gyrus and R. postcentral gyrus 

R. precentral gyrus and R. superior temporal gyrus 

R. superior frontal gyrus and R. superior temporal gyrus 

R. inferior frontal gyrus and R. precentral gyrus 

R. inferior frontal gyrus and R. middle frontal gyrus 

R. precentral gyrus and R. superior frontal gyrus 

R. postcentral gyrus and R. superior frontal gyrus 

L. middle frontal gyrus and R. precentral gyrus 

L. inferior frontal gyrus and R. precentral gyrus 

R. lateral orbitofrontal gyrus and R. precentral gyrus 

L. lateral orbitofrontal gyrus and R. postcentral gyrus 

R. middle orbitofrontal gyrus and R. precentral gyrus 

R. postcentral gyrus and R. superior temporal gyrus 

R. inferior frontal gyrus and R. postcentral gyrus 

L. middle frontal gyrus and R. superior temporal gyrus 

L. middle orbitofrontal gyrus and R. precentral gyrus 

R. middle orbitofrontal gyrus and R. postcentral gyrus 

L. middle frontal gyrus and R. postcentral gyrus 

L. superior frontal gyrus and R. inferior frontal gyrus 

Regions with higher coherence in ASD 

L. angular gyrus and L. middle occipital gyrus 

L. superior occipital gyrus and L. supramarginal gyrus 

L. middle occipital gyrus and L. middle temporal gyrus 

L. angular gyrus and R. middle occipital gyrus 

L. middle occipital gyrus and R. middle temporal gyrus 

L. superior occipital gyrus and R. inferior occipital gyrus 

L. angular gyrus and L. supramarginal gyrus 

L. middle occipital gyrus and L. superior occipital gyrus 

L. middle occipital gyrus and L. superior temporal gyrus 

L. middle temporal gyrus and L. supramarginal gyrus 

(C) Brain region pairs (pathway) for 30 to 45 Hz 

Regions with higher coherence in NT 

L. lateral orbitofrontal gyrus and R. precentral gyrus 

L. middle frontal gyrus and R. precentral gyrus 

R. middle frontal gyrus and R. superior temporal gyrus 

L. inferior temporal gyrus and R. precentral gyrus 



-5.18 

-5.18 

-4.58 

-4.55 

-4.21 

-4.13 

-4.09 

-3.94 

-3.92 

-3.9 

-3.9 

-3.87 

-3.84 

-3.76 

-3.58 

-3.55 

-3.48 

-3.32 

-3.3 

-3.26 

-3.22 

-3.22 

-3.17 

-3.16 

-3.14 

-3.1 

-3.07 

3.66 
3.65 
3.58 
3.51 
3.48 
3.21 
3.19 
3.15 
3.08 
3.08 
z-score 

-4.04 
-3.99 
-3.86 
-3.79 



<0.0001 

0.0025 

<0.0001 

<0.0001 

<0.0001 

<0.0001 

<0.0001 

<0.0001 

<0.0001 

<0.0001 

0.0001 

0.0001 

0.0001 

0.0003 

0.0005 

0.0005 

0.0007 

0.0009 

0.0001 

0.0017 

0.0013 

0.0014 

0.0017 

0.0022 

0.0024 

0.0021 

0.0022 

0.0003 
0.0003 
0.0004 
0.0005 
0.0006 
0.0014 
0.0014 
0.0017 
0.0022 
0.0022 
P value 

0.0001 
0.0001 
0.0002 
0.0002 
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Table 4 Brain region pairs with statistically significant between group differences in coherence for averted gaze in the 
(A) low (0 to 15 Hz), (B) beta (15 to 30 Hz), and (C) gamma (30 to 45 Hz) frequency bands (Continued) 

L. middle frontal gyrus and R. postcentral gyrus -3.76 0.0003 

R. middle frontal gyrus and R. postcentral gyrus -3.74 0.0005 

R. middle frontal gyrus and R. precentral gyrus -3.72 0.0003 

R. inferior frontal gyrus and R. middle frontal gyrus -3.64 0.0003 

L. middle frontal gyrus and R. superior temporal gyrus -3.6 0.0003 

L. inferior frontal gyrus and R. precentral gyrus -3.6 0.0004 

L. middle orbitofrontal gyrus and R. precentral gyrus -3.59 0.0004 

L. lateral orbitofrontal gyrus and R. postcentral gyrus -3.53 0.0006 

L. superior frontal gyrus and R. precentral gyrus -3.46 0.0007 

L. superior frontal gyrus and R. postcentral gyrus -3.46 0.001 

L. inferior frontal gyrus and R. postcentral gyrus -3.44 0.0009 

R. inferior frontal gyrus and R. precentral gyrus -3.4 0.0007 

L. lateral orbitofrontal gyrus and R. inferior frontal gyrus -3.38 0.0007 

R. inferior frontal gyrus and R. superior temporal gyrus -3.37 0.0008 

L. lateral orbitofrontal gyrus and R. middle frontal gyrus -3.32 0.0009 

L. superior frontal gyrus and R. superior temporal gyrus -3.29 0.001 

L. middle frontal gyrus and R. inferior frontal gyrus -3.26 0.001 1 

L. inferior temporal gyrus and R. postcentral gyrus -3.2 0.0014 

R. inferior frontal gyrus and R. postcentral gyrus -3.2 0.0016 

L. middle orbitofrontal gyrus and L. precentral gyrus -3.17 0.0017 

L. middle orbitofrontal gyrus and R. postcentral gyrus -3.16 0.002 

L. superior frontal gyrus and R. inferior frontal gyrus -3.03 0.0024 

L. inferior temporal gyrus and R. inferior frontal gyrus -2.97 0.003 

L. precentral gyrus and R. superior temporal gyrus -2.97 0.0031 

L. inferior frontal gyrus and R. superior temporal gyrus -2.91 0.0046 

L. inferior frontal gyrus and R. inferior frontal gyrus -2.84 0.0047 

L. lateral orbitofrontal gyrus and R. superior temporal gyrus -2.83 0.0047 
Regions with higher coherence in ASD 

R. angular gyrus and R. fusiform gyrus 4.38 <0.0001 

L. angular gyrus and R. angular gyrus 4.13 <0.0001 

L. superior occipital gyrus and R. angular gyrus 3.92 0.0001 

R. angular gyrus and R. cuneus 3.91 0.0001 

R. fusiform gyrus and R. supramarginal gyrus 3.87 0.0002 

L. middle occipital gyrus and R. angular gyrus 3.8 0.0001 

L. superior parietal gyrus and R. fusiform gyrus 3.71 0.0003 

R. angular gyrus and R. inferior temporal gyrus 3.67 0.0002 

L. lingual gyrus and R. angular gyrus 3.61 0.0003 

R. fusiform gyrus and R. middle occipital gyrus 3.59 0.0004 

R. angular gyrus and R. middle occipital gyrus 3.54 0.0004 

R. angular gyrus and R. lingual gyrus 3.49 0.0006 

L. middle occipital gyrus and R. middle occipital gyrus 3.44 0.0007 

L. inferior occipital gyrus and R. angular gyrus 3.39 0.0007 

L. angular gyrus and R. middle occipital gyrus 3.38 0.0007 

L. superior occipital gyrus and R. supramarginal gyrus 3.36 0.0008 
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Table 4 Brain region pairs with statistically significant between group differences in coherence for averted gaze in the 
(A) low (0 to 15 Hz), (B) beta (15 to 30 Hz), and (C) gamma (30 to 45 Hz) frequency bands (Continued) 



R. fusiform gyrus and R. gyrus rectus 


3.34 


0.0014 


L. angular gyrus and L. superior temporal gyrus 


3.32 


0.0009 


L. angular gyrus and R. fusiform gyrus 


3.32 


0.001 


R. cuneus and R. supramarginal gyrus 


3.29 


0.001 


R. angular gyrus and R. inferior occipital gyrus 


3.28 


0.001 


L. angular gyrus and R. supramarginal gyrus 


3.25 


0.0014 


L. superior occipital gyrus and R. superior frontal gyrus 


3.14 


0.0017 


R. inferior temporal gyrus and R. middle occipital gyrus 


3.14 


0.0018 


L. angular gyrus and R. inferior temporal gyrus 


3.12 


0.0018 


R. angular gyrus and R. middle orbitofrontal gyrus 


3.1 


0.002 


L. gyrus rectus and L. superior occipital gyrus 


3.1 


0.0031 


L. angular gyrus and L. lateral orbitofrontal gyrus 


3.06 


0.0022 


L. superior parietal gyrus and R. inferior temporal gyrus 


3.04 


0.0024 


L. superior occipital gyrus and R. middle occipital gyrus 


2.99 


0.0032 


R. angular gyrus and R. lateral orbitofrontal gyrus 


2.98 


0.0028 


L. superior occipital gyrus and R. lateral orbitofrontal gyrus 


2.98 


0.003 


L. middle occipital gyrus and R. fusiform gyrus 


2.93 


0.0034 


R. gyrus rectus and R. inferior temporal gyrus 


2.93 


0.0042 


L. angular gyrus and L. middle occipital gyrus 


2.91 


0.0037 


L. inferior occipital gyrus and R. fusiform gyrus 


2.9 


0.0038 


L. fusiform gyrus and R. angular gyrus 


2.9 


0.0046 


L. angular gyrus and R. lateral orbitofrontal gyrus 


2.85 


0.0044 



Positive z-scores reflect higher coherence in ASD. Negative z-scores reflect higher coherence in NT. 
**P <0.01. 
***P <0.001 . 
****P <0.0001. 



(P <0.001 - right inferior and superior frontal gyri; P <0.01 - 
right inferior frontal and middle frontal), between right 
frontal and right temporal cortices (P <0.01), and unex- 
pectedly, between bilateral occipital regions (P <0.01) was 
related to lower scores and fewer ASD symptoms in the 
ADI-R. Higher symptom reporting was only related to in- 
creased coherence between the right precentral and right 
orbitofrontal regions. 

In the gamma frequency band, higher coherence be- 
tween bilateral frontal regions (particularly right precen- 
tral gyri) and all other cortical regions bilaterally was 
related to a lower number of ASD symptoms as reported 
on the ADI-R. Higher gamma frequency band coherence 
between the right precentral gyrus and left inferior 
frontal (P <0.001), left lateral orbitofrontal (P <0.005), 
left middle frontal (P <0.01), and the left inferior tem- 
poral cortex (P <0.005) was related to lower scores and 
fewer symptoms on the ADI-R. Higher gamma band co- 
herence between right lateral orbitofrontal and left su- 
perior occipital (P <0.005) and the left angular gyrus 
(P <0.005) was also related to fewer ASD symptoms. Fi- 
nally, higher gamma band activity between the right 



temporal gyrus and left superior frontal gyrus (P <0.01) 
was related to lower scores on the ADI-R. 

An opposite and strongly lateralized right hemispheric 
finding was also obtained when examining gamma activ- 
ity. Higher gamma frequency band coherence between 
the right parietal cortex (that is, angular gyrus) and right 
orbitofrontal {P <0.005), temporal (P <0.005 - fusiform; 
P <0.01 - inferior temporal), and occipital regions 
(P < .001 - lingual and cuneus; P <0.005 - inferior occipi- 
tal) as well as between bilateral occipital cortical regions 
and right temporal cortex (P <0.005 - bilateral fusiform; 
P <0.01 - inferior temporal) was related to a greater 
number of ASD symptoms as reported on the ADI-R. 

Discussion 

Processing eye gaze is a vital ability as it provides socially 
relevant information about ones environment, allows us 
to make inferences about the possible intentions of 
others, and is one of the most important aspects of non- 
verbal communication. Although gaze processing deficits 
are a seminal, early, and enduring behavioral deficit in 
ASD, a comprehensive characterization of the neural 
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Figure 4 Regional cortical differences in synchronous beta and gamma activity during averted gaze. Number of region-to-region cortical 
differences within the low (0 to 15 Hz), beta (15 to 30 Hz), and low gamma (30 to 45 Hz) frequency bands in regions of either statistically higher 
coherence in NT or ASD participants. 



60 r 




Intrahemispheric 
Left 

- Intrahemispheric 
Right 

Interhemispheric 



Figure 5 Intra- and inter-hemsispheric cortical differences in synchronous beta and gamma activity during averted gaze. Number of 
intra- and inter-hemispheric cortical differences within the low (0 to 15 Hz), beta (15 to 30 Hz), and low gamma (30 to 45 Hz) frequency bands in 
regions of either statistically higher coherence in NT or ASD participants. 
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Table 5 Significant Kendall Tau correlation coefficients between the ADI-R and brain region pairs with statistically 
significant between group differences in coherence for averted gaze for each frequency band and direct gaze 
collapsed across frequency band 



Averted gaze 












Direct gaze 


Oto 15 Hz 


T 


15 to 30 Hz 


T 


30 to 45 Hz 


T 


Collapsed across all t 
frequencies 


L. angular gyrus and L. 
lingual gyrus 


0.55*** 


L. lateral orbitofrontal 
gyrus and R. postcentral 
gyrus 


-0.53** 


L. angular gyrus and R. 
lateral orbitofrontal gyrus 


-0.55*** 


L. angular gyrus and L. 0.46*** 
supramarginal gyrus 


L. angular gyrus and L. 
middle temporal gyrus 


0.56*** 


L. middle frontal gyrus 
and R precentral gyrus 


-0.60*** 


L. inferior frontal gyrus and 
R. precentral gyrus 


-0.66**** 


L. inferior frontal gyrus -0.41** 
and R. precentral gyrus 


L. angular gyrus and L. 
superior temporal gyrus 


U.DU 


L superior occipital gyrus 
and R. inferior occipital 
gyrus 


-0.52** 


L. inferior occipital gyrus 
and R. fusiform gyrus 


U.j/ 


L. lateral orbitotrontal -U.44 
gyrus and R. precentral 
gyrus 


L. angular gyrus and R. 
angular gyrus 


0.54*** 


R. inferior frontal gyrus 
and R. middle frontal 
gyrus 


-0.51** 


L. inferior temporal gyrus 
and R. precentral gyrus 


-0.59*** 




L. angular gyrus and R. 
middle temporal gyrus 


0.50** 


R. inferior frontal gyrus 
and R. superior frontal 
gyrus 


-0.63**** 


L. lateral orbitofrontal gyrus 
and R. precentral gyrus 


-0.57*** 




L. inferior frontal gyrus 
and R. middle frontal 
gyrus 


-0.56*** 


R. middle orbitofrontal 
gyrus and R. precentral 
gyrus 


0.59*** 


L. lingual gyrus and R. 
angular gyrus 


0.57*** 




L. inferior temporal gyrus 
and R. cuneus 


0.66**** 


R. superior frontal gyrus 
and R. superior temporal 
gyrus 


-0.52** 


L. middle frontal gyrus and 
R. precentral gyrus 


-0.53** 




L. middle occipital gyrus 
and L. middle temporal 
gyrus 


0.52** 






L. middle occipital gyrus 
and R. fusiform gyrus 


0.51** 




R. angular gyrus and R. 
superior frontal gyrus 


-0.56*** 






L. precentral gyrus and R. 
superior temporal gyrus 


-0.50** 




R cuneus and R. gyrus 
rectus 


-0.59*** 






L. superior frontal gyrus 
and R. superior temporal 
gyrus 


-0.50** 




R middle frontal gyrus 








L. superior occipital gyrus 


-0.60*** 





gyrus 



and R. lateral orbitofrontal 
gyrus 

R. angular gyrus and R. 0.64*" 
cuneus 

R. angular gyrus and R. 0.62* : 
fusiform gyrus 

R. angular gyrus and R. 0.51*" 
inferior occipital gyrus 

R. angular gyrus and R. 0.53*" 
inferior temporal gyrus 

R. angular gyrus and R. 0.64*" 
lingual gyrus 

R. angular gyrus and R. 0.71*" 
middle occipital gyrus 

R. angular gyrus and R. 0.56*" 
middle orbitofrontal gyrus 

R. fusiform gyrus and R. 0.61*' 
middle occipital gyrus 

R. inferior temporal gyrus 0.53*" 
and R. middle occipital gyrus 



*P <0.01 . 
**P <0.005. 
***P <0.001 . 
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processes mediating abnormal gaze processing in ASD 
has yet to be conducted. 

MEG studies exploring other well-described behavioral 
phenomena, such as deficits in face and emotion process- 
ing [47-50] , have recently been reported following the 
growth of more unified theories that include components 
of altered connectivity, an imbalance in excitatory to 
inhibitory neural transmission, and impaired neural syn- 
chrony as fundamental pathophysiological mechanisms of 
ASD. These preliminary studies have reported findings 
suggestive of abnormal functional organization, aberrant 
pathway development, and possible altered hemispheric 
specialization. 

Recent investigations have also explored neural synchro- 
nization (phase coherence) and connectivity during both 
face [33,51] and auditory/language processing in ASD in 
an attempt to better elucidate aberrant patterns of con- 
nectivity and identify potential biomarkers [52-56], with a 
specific focus on gamma power and oscillatory activity. 
The results are generally equivocal with studies reporting 
higher regionally induced gamma power, higher and lower 
regionally evoked gamma power, and reductions in phase 
consistency or timing of power (that is, phase locking). 
This may be partially accounted for by the rather signifi- 
cant heterogeneity in this population coupled with the 
extreme variability in the methodologies being used and 
populations being examined. 

Resting state investigations exploring oscillatory activity 
have only very recently emerged [57-59] consistent with 
efforts to more fully characterize global aberrant connect- 
ivity patterns in ASD. Indeed, Tsiaras et al [59], reported 
finding attenuated short-range connectivity in adults with 
ASD within bilateral temporal and frontal regions and left 
parietal regions, although significant differences between 
specific frequency bands were not apparent. In a study of 
resting state in children with ASD compared to typically 
developing children, Cornew and colleagues [57] recently 
reported finding various differences in oscillatory activity 
including increased theta and alpha power in parietal and 
occipital regions and additionally increased alpha power 
in temporal regions in ASD. The authors further reported 
finding greater relative delta in right frontal regions in ASD. 
These results are generally consistent with neurophysio- 
logical findings from EEG [60], but they fail to bring us 
closer to a comprehensive understanding the relationship 
between the behavioral phenotype and the pathophysio- 
logical mechanisms of this disorder. 

Given its central phenotypic prominence [19-23] in 
ASD, gaze processing clearly surfaces as a strong and 
enduring endophenotypic candidate. In order to more 
broadly understand the relationship between a core and 
enduring behavioral deficit in ASD and its neuropathol- 
ogy, we hypothesized a priori that processing eye gaze in- 
formation was likely to more precisely characterize 



aberrant beta and gamma band oscillatory activity and 
potentially aberrant connectivity. Our results are very con- 
sistent with the aforementioned investigations and 
revealed that ASD participants demonstrated lower coher- 
ence between bilateral frontal (particularly right frontal) 
and right pre-and postcentral regions and superior tem- 
poral regions when passively viewing gaze. In contrast, 
ASD participants demonstrated higher coherence between 
sensory association cortices (that is, temporo-parieto- 
occipital) in all frequency bands, particularly within the 
low frequency range, as well as in those associated with 
both short- and long-range transmission. 

Our results also reveal a very clear relationship between 
aberrant oscillatory activity and elevated ASD symptoms, 
specifically an increase in low frequency activity posterio- 
rally and bilaterally as well as an increase in gamma activ- 
ity in right posterior temporo-parietal-occipital regions. In 
contrast, lower symptomatology appears related to 
increased low frequency coherent activity between the 
right frontal cortex with left frontal and right parietal and 
occipital regions as well as related to increased oscillatory 
activity in the gamma frequency band between right 
frontal regions with all other left hemisphere lobes 
(frontal, temporal, parietal, and occipital) and between the 
left frontal and right temporal lobes (see Figure 6). The 
regional brain differences noted in ASD participants 
displaying fewer clinical features are very consistent with 
the purported brain regions known to underlie social 
cognition. 

Disruptions in neuroanatomical pathways 

Recent multi-modal imaging methods have started to 
illuminate networks involved in direct and averted gaze 
in NT using, for example, combined fMRI-DTI [61]. Dy- 
namic gaze shifts have been found to increase activation 
in a well-established region of the social network includ- 
ing the right pSTS, anterior insula, and fusiform gyrus 
with direct connections noted between the right pSTS 
and anterior insula. These regions are thought to be sup- 
ported by long-range network connections that project 
via the superior longitudinal fasciculus and are believed 
to be critical for extracting social meaning of eye gaze 
shifts [61]. Consistent with this purported social net- 
work, compared to those with ASD, our NT partici- 
pants demonstrated significantly higher coherence 
between the right middle and inferior frontal gyri and 
the right superior temporal regions within the beta 
band frequency known to be critical for long-range 
connectivity. 

In contrast, our ASD participants demonstrated higher 
gamma power in inter-hemispheric connections between 
the left and right parietal lobes and intra-hemispherically 
between the right parietal lobe (angular gyrus) and tem- 
poral regions. The angular gyrus has been implicated in 
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Figure 6 Relationships between oscillatory activity in the low, beta, and low gamma frequency bands during averted gaze and autism 
symptomatology measured by the ADI-R. Note increased low and gamma frequency band oscillatory connections in posterior temporo-parieto- 
occipital regions associated with a higher number of ASD symptoms while high frequency gamma oscillatory connections between right frontal 
regions and all left hemisphere lobes is associated with fewer ASD symptoms. 



a number of processes including reading and number 
comprehension, numerical processing, visual attention, 
and social cognition. It is a cross-modal region where sen- 
sory information from the visual, auditory, and tactile 
senses converge allowing for a combined and integrated 
percept. It is essential for the manipulation of mental rep- 
resentations and reorienting of attention [62]. We believe 
that heightened connectivity in these regions without 
appropriate regulatory or contextual feedback from 
frontal regions may result in exquisite sensory sensitivity, 



acceleration of letter, number, and word recognition with 
limited comprehension or applied skills, or an over- 
allocation of attention to information without a clear ap- 
preciation of its relevance; a neurocognitive pattern often 
noted in ASD. A strengthening of connectivity between 
the right angular gyrus and inferior temporal regions with- 
out frontal mediation, particularly from medial prefrontal 
regions, may further contribute to a heightened attention 
to the invariant features of the face or its components 
without an ability to extract essential social relevance. Our 
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findings are consistent with models of ASD proposed by 
Brock and colleagues [25] as well as Belmonte et al [27] 
who have suggested that abnormally high levels of high- 
frequency neural activity and over connectivity within lo- 
calized brain regions causes impaired discrimination of 
brain signal from background brain noise. Interestingly 
and unexpectedly, our ASD participants demonstrated a 
greater number of low frequency cross-hemispheric con- 
nections, and particularly in posterior cortical regions. The 
results suggest that conceptualizing ASD as a disorder of 
either over- or under-connectivity in long- or short-range 
connections may not fully capture the pathophysiology. 

Synchronous oscillatory activity during gaze cueing as a 
biomarker of ASD? 

Approximately 200 investigations have been published 
proposing neuroanatomical markers of ASD, although 
the results have often been in conflict or unreplicated 
[63]. The majority of this work has been in older children 
and adults using structural imaging. Recently, Ecker and 
colleagues [64] used linear support vector machines 
(SVM), a machine learning method that identifies patterns 
in data by identifying hyperplanes that maximally differen- 
tiate categories or groups, classified 20 adults with and 
without ASD. They used five structural classifiers of cor- 
tical gray matter to correctly classify 85% of those with 
ASD, with 90% sensitivity and 85% specificity. Left hemi- 
spheric classifiers were more accurate than right, and 
cortical thickness was the best classifier. Hemisphere 
laterality is an area that remains relatively unexplored in 
ASD, and with their reported method the authors were 
unable to determine if individuals with ASD displayed a 
higher (lower) degree of cortical asymmetry. Similarly, a 
recent fMRI investigation by Dinstein and colleagues [29] 
of sleeping toddlers with ASD revealed significantly 
weaker inter-hemispheric synchronization (that is, weak 
'functional connectivity across the two hemispheres) in in- 
ferior frontal (IFG) and superior temporal (STG) regions 
for which early over-lateralization of language function 
was suggested. However, again, directionality of lateraliza- 
tion to the left or right hemisphere could not be deter- 
mined from their data. Our method allows us to provide a 
direct numerical comparison between pathways, both 
inter- and intra-hemispherically, and to examine group 
differences between frequency bands known to underlie 
short- and long-range connectivity. Given the passive na- 
ture of the task, this methodology could be easily applied 
to preschoolers. Only 18% of children with ASD are iden- 
tified by the age of 3 years; even later for children with 
milder forms (average, 6.3 years) [65]; well beyond when 
children can benefit maximally from early intervention. 
This method allowed coherence to be imaged in source 
rather than sensor space, achieving a significant degree of 
disentanglement of signals. Moreover, it provided better 



resolution and precision of the underlying networks gen- 
erating the signals since coherence in sensor space is 
smeared due to current spread as described by Srinivasan 
et al [66]. Neural synchrony between the frontal cortex 
and pre- and postcentral gyrus, between frontal and 
superior temporal cortex in the beta band, and between 
angular and fusiform gyri in the gamma band were 
particularly discrepant between the groups. 

We recognize that there are a number of weaknesses 
with this pilot study. First, the sample size of this pilot 
phase was small which limits the external validity and 
conclusions that can be drawn. Although the findings 
were quite robust, a replication of this study with a larger 
number of participants is necessary. Moreover, the results 
suggest a strong need for funding of large scale and longi- 
tudinal studies with this technology in ASD. 

Second, our ASD participants remained on their psy- 
chotropic medication regimens during the study given the 
perceived cost-benefit of motion artifact for participants 
taking ADHD medications as well as the practical chal- 
lenges associated with withdrawing from medications with 
longer half-lives. However, we recognize that when one is 
evaluating a possible pathophysiological mechanism in 
any psychiatric or neurodevelopmental disorder, a vital 
issue is whether the observed neural abnormality is 
present at the outset of the disorder and predates medica- 
tion exposure and the deleterious effects of long-term 
illness. As noted, five of our participants were on one or 
more of the following medications: psychostimulants, 
antidpressants, anxiolytics, antipsychotics, or mood stabi- 
lizers. Recent investigations exploring alterations in oscil- 
latory activity in populations taking medications similar to 
our cohort (for example, ADHD, bipolar disorder, and 
schizophrenia) have reported variable results with respect 
to medication status on oscillatory activity and the 
findings to date have been primarily reported in adult 
populations [67-69]. Wilson and colleagues [69] recently 
examined time estimation and oscillatory activity using 
MEG in medicated and unmedicated individuals with 
ADHD. Relative to controls, unmedicated participants 
were reported to display less accurate time estimation and 
weaker gamma activity in frontal cortex, specifically anter- 
ior cingulate, supplemental motor areas, and the left pre- 
frontal cortex. Following medication administration, the 
patients demonstrated small but significant increases in 
gamma-band activity across the same neural regions, 
which was related to improved time estimation accuracy. 
Exploratory analysis also revealed stronger delta activity in 
frontal regions in the control participants relative to those 
with ADHD, regardless of medication status, possibly sug- 
gesting that such alterations in delta activity may not be re- 
sponsive to stimulant medications. Although we recognize 
it is difficult to compare across populations, recent investi- 
gations have demonstrated high co-morbidity with ADHD 
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and ASD, suggesting similar aberrant neural substrates. 
We further reviewed the literature to explore possible 
effects of antipsychotics on oscillatory power. Minzenberg 
and colleagues [67] recently reported reductions in gamma 
oscillatory power during an executive task that was inde- 
pendent of medication status in first-episode schizophre- 
nia. Fifty-three first episode patients with schizophrenia 
(21 without antipsychotic treatment), aged 13 to 30 years, 
underwent EEG while performing a cognitive control task. 
Both the medicated and unmedicated patient subgroups 
were impaired on the behavioral task and displayed lower 
frontal induced gamma power. In contrast, there were no 
significant group differences in theta power between con- 
trols or the medicated or unmedicated patient subgroups. 
Given the above, we might have expected the effect of psy- 
chostimulants and antipsychotics to attenuate the frontal 
gamma frequency band group differences noted between 
our groups. Ozerdem et al [68] examined alpha and beta 
oscillatory activity during a visual odd-ball paradigm and 
EEG in 10 bipolar patients before and after valproate treat- 
ment. At baseline, drug-free individuals with bipolar 
disorder in the hypomanic or manic phase of illness dem- 
onstrated aberrant alpha and beta oscillatory activity to vis- 
ual target stimuli compared to healthy controls. Individuals 
with bipolar disorder displayed significantly higher beta 
activity in occipital cortices compared to controls, and 
patients were devoid of an occipital-frontal alpha domin- 
ance that was noted in the control group. Following treat- 
ment with valproate, occipital beta responses dropped to 
control levels, and a further and significant decrease in 
occipital alpha and unchanged frontal alpha response was 
noted in patients with bipolar disorder. The authors 
suggest that the latter may suggest the inhibitory effect of 
valproate via GABAergic inhibition. That is, inhibition of 
occipital alpha activity while trying to manage cortical 
hyperactivity may be an unwanted effect of valproate use. 
Similarly, if the effect of valproate is an attenuation of oc- 
cipital alpha activity, the impact of this medication would 
have also been to decrease our current effect of increased 
low frequency band activity in posterior regions. However, 
these assumptions should be interpreted cautiously as it is 
important to remember that investigations examining 
auditory and visual evoked oscillatory activity in humans 
have shown that the evoked responses are topography and 
stimulus modality dependent and that brain structures 
with different resonance properties may be dependent on 
the stimuli. 

Finally, with respect to the frequency ranges that were 
examined, while we did not examine activity greater 
than 45 Hz, we recognize that the gamma band fre- 
quency extends to approximately 100 Hz. For this pilot 
project, we chose to remain below 60 Hz given that 
group differences in gamma power are the most robust 
in the 40 Hz range and in order to avoid additional 



artifact. Despite these limitations, this pilot study provided 
the first cogent understanding of whole-brain patterns of 
coherence that underlie direct and averted eye gaze in 
ASD. 

Conclusions 

One of the hallmarks of ASD is a failure to detect and/ 
or respond in a typical manner to information conveyed 
by eye gaze. This study characterized whole-brain pat- 
terns of synchrony in ASD compared to NT during 
direct and averted eye gaze processing while undergoing 
MEG. Results revealed: (1) higher coherence and 
synchronization in temporo-parietal-occipital brain 
regions across all frequencies in ASD, particularly within 
the low frequency range; (2) a higher number of low 
frequency cross-hemispheric coherent connections; and 
(3) a near absence of right intra-hemispheric synchrony 
in the beta frequency band in ASD. This preliminary 
examination of the relationship between neural syn- 
chrony and ASD symptomatology further revealed a 
unique pattern of relatively mutually exclusive findings: 
(1) higher synchronization in the low frequency band (0 
to 15 Hz) between left temporo-parieto-occipital regions 
and higher synchronization in the low gamma frequency 
band (30 to 45 Hz) in right temporo-parieto-occipital 
regions was associated with more severe ASD symptom- 
atology while (2) higher synchronization in the low 
frequency band between right front-parieto-occiptal 
regions as well as in the gamma frequency band between 
bilateral frontal and right frontal with left temporo- 
parieto-occipito regions was associated with fewer ASD 
symptoms. This altered pattern of oscillatory activity 
may contribute to aberrant connectivity that underlies 
the failure of individuals with ASD to appropriately 
orient to eye gaze, which has a cascading negative effect 
on typical social and language development. 
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